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,i__t ■ Abstract. We derive the asymptotic maximum- likelihood phase estimation 

^ , uncertainty for any interferometric protocol where the positions of the probe particles 

^ 1 are measured to infer the phase, but where correlations between the particles are not 

O" 1 accessible. First, we apply our formula to the estimation of the phase acquired in the 

Mach-Zehnder interferometer and recover the well-know momentum formula for the 
04 1 phase sensitivity. Then, we apply our results to interferometers with two spatially 

separated modes, which could be implemented with a Bose-Einstein condensate 
trapped in a double-well potential. We show that in a simple protocol which estimates 
the phase from an interference pattern a sub shot-noise phase uncertainty of up to 
can be achieved. One important property of this estimation protocol 
00 1 is that its sensitivity does not depend on the value of the phase 8, contrary to the 

sensitivity given by the momentum formula for the Mach-Zehnder transformation. 
Finally, we study the experimental implementation of the above protocol in detail, by 
• • ■ numerically simulating the full statistics as well as by considering the main sources of 

. t-H . detection noise, and argue that the shot-noise limit could be surpassed with current 

^ " technology. 
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1. Introduction 

Quantum interferometry, a rapidly developing branch of modern physics, exploits some 
striking features of quantum mechanics in order to build ultraprecise measuring devices 
pp. It employs non-classical states of light EJ HJ [5] or matter EJ EJ El EEl [EH E2] 
to beat the shot-noise limit (SNL) A9 = 1/ \/~N - the limit of precision of parameter 
estimation of a phase shift 9 set by classical physics. Here N is the number of particles 
in the probe state. Recent achievements in the preparation of entangled states of atoms 
have put atomic interferometry [13] in the avant-garde for precise determination of 
electromagnetic [HI [151 05] and gravitational [T71 HH [19] interactions. 

A paradigmatic example and benchmark for every interferometer is the Mach- 
Zehnder interferometer (MZI), where two "beams" propagating along separated paths 
accumulate a relative phase 9 (to be estimated) and are subsequently recombined 
through a beam splitter. A similar protocol which uses internal states instead of 
separate paths is known as Ramsey-spectroscopy [20] • It has been recently realized 
with a Bose-Einstein condensate (BEC), employing two hyperfme states of the atoms 
[9] . The beam-splitter was mimicked by coupling the two modes with a micro- wave pulse 
for a precisely chosen amount of time. Atomic interferometers with spatially separated 
modes, which could for instance be used to measure forces decaying with the distance are 
more challenging in implementation due to the difficulty of performing the beam-splitter 
transformation J2TJ [22] . 

In this paper we consider a simple "double-slit" interferometer in which two beams 
are recombined through a free expansion of two initially localized clouds with the output 
signal obtained by measuring the positions of the particles. After preparation of a 
suitably entangled input state, sub shot-noise (SSN) interferometry requires thus refined 
particle detection at the output. Some new techniques of atom-position measurements, 
like the micro-channel plate [23], the tapered fiber [24] . the light-sheet method [25] or 
techniques involving atomic fluorescence from the lattice [26], give hope for an almost 
100% efficient single-atom detection in the nearby future. Such a tool in principle could 
even give access to atom-atom correlations at all orders. 

It has been shown [271 I2H] that the measurement of the iV-th order correlation 
function is the best possible estimation strategy for inferring the phase between 
two interfering BEC wave-packets, and allows to reach the Heisenberg limit of 
phase uncertainty. However, even with small BECs, the measurement of the iV-th 
order correlation function would require substantial experimental effort, since a huge 
configurational space must be probed with sufficient signal to noise ratio. 

The difficulty of measuring high order correlation functions is the motivation for this 
work. We show that by measuring positions of particles the phase can be estimated with 
a SSN phase uncertainty using the simple one-body density. We consider two possible 
detection scenarios: the output signal consists of (a) the positions of single atoms, or 
(b) the number of atoms per pixel, which corresponds to the commonly employed least- 
squares estimation from the fit to the average density [29] . For both cases, we compute 
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Figure 1. Schematic representation of the atom position measurement, where the 
detectors (boxes) turn yellow (no. 1, 3 and 7) when they are hit by atoms (spheres). 
Some detectors however remain gray (as the detector no. 8), so the detection efficiency 
drops, and sometimes a neighbouring box turns yellow (no. 5 instead of no. 4), which 
limits the spatial resolution. 

the general asymptotic phase uncertainty of a maximum-likelihood phase estimation 
scheme using the one-body density only. 

We verify that in an MZI situation, estimating the phase from the one-body density 
is equivalent to estimation from the average population imbalance between the arms of 
the interferometer, and recover the known result that the sensitivity can saturate the 
Heisenberg limit. 

We then apply this estimation protocol to our case of interest, namely, two 
interfering BEC wavepackets. We provide an analytical expression for the phase 
uncertainty which shows SSN scaling with phase squeezed states in input. Contrary 
to the MZI, the phase uncertainty does not depend on the phase 9, which could be an 
important advantage. We then analyze the full statistics of the phase estimation by 
numerically simulating experiments with a realistic number of particles. We find that 
already small statistical samples are sufficient to saturate the analytical asymptotic 
prediction for the phase uncertainty. Finally, we discuss the main sources of noise 
affecting the interferometric precision, which we expect to concern the atom detection 
stage. We argue that including the effect of imperfect detection the SNL can still be 
surpassed, and that the necessary amount of squeezing could be achieved in a realistic 
double-well setup. 

The paper is organized as follows. In Section [2] we derive the general expression 
for the uncertainty of the phase estimated from the one-body density. We consider 
two possible detection scenarios. We assume either having access to positions of single 
particles or a coarse grained measurement due to a limited detection capability. In 
Section EJ we apply the above results to the case of the MZI. In Section @] we turn to 
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our case study where the interferometer consists of a simple phase-imprint between the 
modes followed by the ballisitic expansion of the mode functions. We also show how the 
phase uncertainty for the phase estimation from the density of the inteference pattern 
is influenced by characteristic sources of noise. Some details of analytical calculations 
are presented in the Appendices. 



2. Sensitivity of the one-body density estimator 

An interferometric estimation protocol can be divided in: (i) the interferometer 
transformation, (ii) the measurement at the output, and (iii) the phase inference through 
an estimator. In (i), the interferometer imprints a relative phase 9 on the iV-particle 
input state l^in)- Such transformation can be represented by a unitary evolution 
operator U(9) = e l6h , where h is linear - i.e. can be written as a sum of operators acting 
on each particle separately - and does not depend on 9. In the Heisenberg picture, 
the field operator evolves as &(x\9) = U\9)^ (x)U (9) and |^>; n ) is unchanged. In the 
following, we will specify our arguments to the case of particle position measurements 
at the output, although the results of Sec. I2.1l are valid in general. We will thus consider 
next that in stage (ii), upon leaving the interferometer, the positions of the particles are 
detected. As observed above, due to the difficulty of obtaining high order correlation 
functions, often only the lowest, namely the density, can be precisely measured. This 
density, when normalized, gives 



the probability density of measuring a single particle at position x given 9. The average 
value is calculated for the input state iV'm)- Next, we assume having no access to the 
correlations, so the phase is inferred only using pi(x\9) and the measurement outcomes 
obtained in stage (ii). 



2.1. Single-atom detection 

When in (ii) single-atom detection is performed, see Fig. HJ each experiment gives the 
positions of the N atoms, x^ = x\ , . . . ,x$, where the label % = 1, . . . ,m indicates 
the particular iteration of the experiment. The accumulated set of data, iV times m 
positions, is used to construct the likelihood function, 

m N 

^)=nn^V). (2) 

1=1 k=l 

The prescription of the Maximum Likelihood Estimator (MLE) is to infer the phase 
9jj![£ in stage (iii) as the value of ip which maximizes Eq. (J2]). As demonstrated in the 



Appendix A[ this MLE is consistent, i.e 9y^ — > 9 for m — > oo. There, we also obtain 



the uncertainty of the estimator equal to 



^---^(l + ^-l)^]. (3) 
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where C depends on the two-body probability density of detecting one particle at x\ 
and the other at x 2 , 

P2{X 1 ,x 2 \o) — iV(iV^T) ' 

and reads [30] 



C= f dx ± dx 2 Vi (gi ,x 2 1 0) dePl ^ ^j 2 ]^ . (5) 
Furthermore Fi is the Fisher information calculated with the probability pi , 

Fi= LJdeM^ 

Equation (j3J) is the first important result of this paper. The phase, which is estimated 
from the single-body probability, depends on both p\ and p 2 . If C = or neglected as in 
[27], Eq. ([3]) provides a shot-noise limited phase uncertainty, since Fi < 1 [27J. Yet, C 
can assume negative values, allowing for SSN phase uncertainty, as will be demonstrated 
below. 

We also underline the generality of the above result, which is valid for any quantum 
state, where the parameter 9 is estimated from the one-body density. In analogy, if two- 
body correlations can be measured in an experiment, four-body correlations would enter 
in the corresponding expression for the asymptotic phase uncertainty. 

2.2. Multiple-atom detection 

It is also important to consider the possibility that the detectors cannot resolve positions 
of individual particles. In this case, we need to use the coarse-grained density and 
assume that in stage (ii) in the z-th experiment the number of atoms n£ in each of the 
k — 1 . . . ribin bins is measured (the bin size Ax must be small to precisely sample the 
density variations). The measurement is repeated m times and the phase is estimated 
from a least square fit of the one-body density to the accumulated data. As discussed 
in detail in [29], the phase uncertainty of such fit is equivalent to the phase uncertainty 
of the MLE with the likelihood function C^{ip) = Y[k=i P(^k\v) , where p(n~k\(p) is the 
probability for detecting — — Yl^i n t^ a t° ms on average in the k-ih bin. For large 
m, according to the Central Limit Theorem, the probability p(nk\f) is a Gaussian with 
a mean (nfc) = AxNpi(xk\<f) and Poissonian fluctuations, A 2 n& = (rife) [27J. For the 
phase uncertainty of this MLE, we obtain similarly as in Section 12.11 that 

E"bin [de{n k )] 2 , v^ n bin n 2 dg{n k ) d B {ni) 
w , fc=l (n k ) + L^=1°*!,I (n k ) (m) ,_v 

m 



2 

k=l (n k ) 



^"bin [dg(n h )] 



where 



alj = (Ax) 2 N [{N - l)p 2 {x k ,xi\6) - Np 1 {x k \e)p 1 {x l \i 
In the continuous limit Ax — > 0, the formulas and flTJ) coincide. 
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3. Estimation from the one-body density with the MZI 

As a benchmark for the phase estimation protocol introduced in Sec. EJ we first consider 
the case where in stage (i) the system acquires the phase 9 in a MZI. The interferometric 
sequence of the MZI consists of three steps. First, the initial two-mode state |^ n ) passes 
a beam-splitter. Then, a relative phase 9 is imprinted between the modes. In the final 
stage of the interferometer, another beam-splitter acts on the state. 

It can be easily shown that the evolution operator of the whole MZI sequence is 
U{9) = e~ %ejy [31]. It is convenient to switch to the Heisenberg picture. The initial field 
operator reads ^>{x) = ip a (x)a + i/jb(x)b, where a/b annihilates a particle from mode a/b 
and the corresponding spatial mode function is i/j a /b(%)- When passing the MZI, this 
operator is transformed as follows 



#(0) = e ieJ ^(x)e- i9j y = ijj a (x) 



9 i ■ 6 
a cos o sin - 

2 2 



i 9^.9 
b cos - + a sin — 



(8) 



We use this result to calculate the one-body [cf. Eq. ([I])] and two-body [cf. Eq. (@|] 
probability densities that enter Eq. (J3J), assuming that both the mode functions are 
point-like and trapped in separate arms of the interferometer, i.e. ip a {x)ipi,(x) = for 
all x. Using an initial iV-particle two-mode state \ipm) = J2j c j\h N — j), where \cj\ 2 is 
the probability of having j atoms in mode a and (N — j) in b, we get 

A 2 e (m) = 1 A 2 J z cos 2 9 + A 2 J x sm 2 9 
ML m (4)2 cos 2 



The details of the derivation are presented in |Appendix B[ The above expression 
coincides with the phase uncertainty for the estimation of the phase 9 from the 
average population imbalance between the two arms of the MZI. This coincidence 
can be explained as follows. When atoms are trapped in two separate arms of the 
interferometer, their positions can be directly translated to the number of particles in 
each arm, without any loss of information. It is then not surprising that the estimators 
basing on the average population imbalance and the average density are equivalent. 
Note that for a particular case of 9 = 0, we get 

= (10) 



where £ n = x/Ny^h is the spin-squeezing parameter [201 E2] related to number 
squeezing of the initial state. The expression (ITUj) can provide up to the Heisenberg 
scaling of the phase uncertainty, once the interferometer is fed with a strongly number- 
squeezed state 



4. Estimation from the one-body density with the interference pattern 



We now turn the attention to our case of interest, to which we devote the rest of the 
manuscript. We consider the whole phase acquisition sequence as consiststing of two 
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steps. First is the phase-imprint performed in absence of two-body interactions, which 
in the two- mode picture is represented by the unitary operator e~ l9Jz , and gives 

4f(x\9) = Mx)e^a + M^e'^b. (11) 

Then, the trap is opened. Since the two-body collisions are assumed to be not present, 
the two wave-packets freely spread and interfere. Atom interactions can strongly 
influence the expansion at an early stage, when the density of the clouds is high [SU ES] • 
We assume that initially i> a /b{ x ) have identical shape, but are centered around ±x . 
After long expansion time (in the so-called "far field") ip a /b{ x ) — e T%h ^ ■ r/> (jy), where 

a = (fi is the atomic mass and t is the expansion time) and k — 2p-, -0 is the Fourier 
transform of the wave-packets at t = 0, the same for ip a /b [27] [28] . Note that we have 

■ x 2 

dropped the common factor e 2s 3 ". Some aspect of this simple interferometric sequence, 
which does not require the implementation of a beam splitter, have been discussed in 
[28]. 

The field operator ^/(x\9) gives pi(xi\9) and P2(%i, X2\9) presented in the 



Appendix C[ which are put into Eq. (|3j). The integrals are performed analytically 
assuming that the interference pattern consists of many fringes, giving 

1 



a2/]H _ 



vT 



v 



V 2 



(12) 



where G = JjV^f is the spin-squeezing parameter [3U] related to phase-squeezing 

of the initial state. Also, we have introduced v = jr(J x ) - i-e. the visibility of the 
interference fringes, see Eq. (1C.1I) and below for details. One important property of 
Eq. (IT2"]) is its independece of the actual value of the phase 9. This can be understood 
as follows. The probabilities ([1]) and (j3J) depend on 9 in the same manner, i.e. via a 
sine or cosine function, see Eqs (10. 1 1) and ( lC.2p . Since the integration in (J5J) and 
runs over the whole space, the shift of the trigonometric functions by a common factor 
9 is irrelevant, since within the envelope ip f^-J there are many interference fringes. 
On the other hand, the sensitivity of the MZI ([9]) depends on 9 and has some optimal 
"working points" because the wave-packets ip a /b(x) do not add up to a wide envelope 
and therefore the argument for ^-independence valid for the interference pattern does 
not apply here. Therefore, the ^-independence of Eq. (fT2"l) might be an advantage of 
this estimation protocol with respect to the momentum formula of the MZI, valid when 
the phase is estimated from the one-body density. 

Note that the phase uncertainty of the MZI [Eq. (ITU]) ] would closely resemble 
the above result if the second term of Eq. ( IT2T) were absent. In the former case, 
^ 2 ^ml benefits from the number-squeezing because the estimator is equivalent to 
the average population imbalance between the arms of the interferometer, therefore a 
state with reduced population imbalance fluctuations decreases the uncertainty A 2 6^ . 
Analogously, in the latter case, the decreasing fluctuations of the relative phase between 
the two modes would improve the estimation precision. 



Sub shot-noise interferometry from measurements of the one-body density 



8 



However, the situation complicates due to the presence of the second term in 
Eq. (|12p . Namely, when £^ drops, so does the fringe visibility v and so the amount 
of information about the phase 9 contained in the one-body density degrades. As a 
consequence, the phase uncertainty of the estimation from pi(x\9) declines. 

It is now important to check whether, due to this interplay between improvement 
from the phase-squeezing and deterioration from the loss of visibility, Eq. ffl2|) can give 
SSN phase uncertainty at all. To this end, we calculate A 2 6^ with phase- squeezed 
states for iV = 100, which we generate by computing the ground state of the two- 
mode Hamiltonian H = —J x + ^jJ^ for negative values of the interaction-to-tunelling 
ratio jfj. For a detailed study of the preparation of phase-squeezed states, including 



experimental sources of noise, see [36J. For every value -^rj, we find \ijj- m ) and calculate 
and v. These values are substituted into Eq. (TI2]) and the resulting A 2 #{^ is plotted 
in Fig. [2] function of £a. 

Clearly, the uncertainty drops below the SNL. We notice the presence of an optimal 



0.012 




0.002 Y 



Figure 2. The phase uncertainty (|T2|) - black solid line - as a function of £^ for 
N = 100 particles. The dashed red line is the phase uncertainty from the squeezing 
parameter (in absence of the second term in Eq. (TT2|)). while the dotted blue line is given 
by the inverse of the QFI. The horizontal dashed line denotes the SNL and the vertical 
dashed line indicates the position of the optimal point. The open green circles are the 
result of a numerical experiment (see text) for four values of ^ = 0.44, 0.59, 0.72 and 



0.86 with m = 10, n r 



4000. 



point, where the gain from the spin-squeezing is balanced by the loss of visibility. The 
figure also shows that the phase uncertainty ffl~2"l) does not saturate the bound set by 
the Quantum Fisher Information (QFI), A 2 #qfi = = ^ 4A 2j 

To support these analytical results for the asymptotics, we also study the full 
statistics of the protocol by simulating a phase-estimation experiment with N = 100 
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particles. We generate the input state with a desired amount of phase-squeezing, 
evaluate the full iV-body probability Pn(xn\@), with which we draw a single realization 
yielding the N positions. We repeat the experiment m times and obtain one value of 
the phase 9^ using MLE with Eq. (j2J). This cycle is performed n rep times and the 
variance of the estimator is calculated on the resulting ensemble. The results, plotted 
(empty circles) in Fig. [2] for four values of m = 10, and n rep = 4000, agree with 
the theoretical value calculated with Eq. f fl2|) . An important information in view of 
an experimental implementation is that though formally the MLE saturates ( fl2l) when 
m — > oo, in practice m = 10 is sufficient to reach the bound, as shown in the upper panel 
of Fig. [31 The lower panel of Fig. [3] shows instead the average value of the estimated 
phase plus uncertainties as a function of m. Here we have chosen the true value of 9 to 
be zero, thus the figure show that the estimator is unbiased. 




Figure 3. The upper panel shows that the phase uncertainty of the "experiment" 
(black points) at the optimal point £^ = 0.44 and N = 100 beats the SNL (which 
is equal to mA 2 9 t ^ = 0.01 here, cf. Fig. [2]) already for m = 1 and converges to 
the theoretical value (dashed red line) for m = 10. The lower panel depicts the 
corresponding average plus error bars indicating uncertainties, while the true value of 
6 is indicated by the dashed red line. 



Our next step is to find the best scaling of Eq. ( 1121) with N. In order to get an 



analytical estimate, we model |^i n ) with a Gaussian, see |Appendix D[ and find that at 
the optimal point mA 2 ^ L ' ,opt = 2iV~ 4 / 3 . This prediction is compared with numerical 
results, where for every iV we evaluate the phase uncertainty ( fl2l) at the optimal state. 
As shown in the inset of the Fig. HI the agreement between the numerics and the 
Gaussian approximation is very good. Also, we numerically obtain the scaling of the 
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QFI at the optimal state mA^Qpj = A _4//3 , which differs from A 2 6^' opt just by a 
factor of 2. Next, we discuss some characteristic experimental imperfections which can 
spoil the phase uncertainty ( IT2|) . 



4-1. Impact of detection imperfections 

Single-particle detection, which is the basis of the estimation scheme discussed in 
Sec. I2.1[ is affected by two dominant sources of noise: limited efficiency and finite spatial 
resolution, see Fig. [TJ The former is incorporated by letting the index k run from 1 to 
n < N in Eq. (12"1). In effect, n replaces N in Eq. (J3]). The finite spatial resolution modifies 
instead both C and F\. We implement it by convoluting the probabilities pi(xi\8) and 
P2(xi,x 2 \0) with p(xi\xi) = -7g=g e~( xi ~ x '^ ^ 2|7 ^ (and analogously for x 2 ) - a Gaussian 
probability of detecting an atom at X\ given its true position x[. The convolutions are 
calculated analytically and the phase uncertainty becomes 



A 2 e [m) - — 



ML " mN 



VI - v 2 e- K ' 2 ° 2 + 1 



77 



77 z/ 2 



(13) 



where k = 2|§ and a = y ^ were defined below Eq. ( TTTT) . Above, we assumed that 
A, n ^> 1 , 77 = jj and the tilde denotes the phase uncertainty in presence of errors. 
Note that for 77 = 100% and o = we recover Eq. ( Tl2l) . In Fig. H we plot Eq. (JBJ) 
as a function of a for iV = 100 and 77 = 100%, 90%, 80% and 70%. For instance, with 
resolution of 1/30-th of a fringe (a = 0.033 x ^ ~ ^) and 77 = 90%, the SSN scaling 
is mA 2 ^ L ' ,opt = 1.48iV~ L16 , shown in the inset. Also note that even in presence of the 
noise, the phase uncertainty ([TBI does not depend on 6. 

The two main sources of detection noise affecting the least-squares fit estimation 
protocol defined in Sec. 12.21 not shown in Fig. [Tj, are imperfect atom counting 
and the finite bin size. To model the former, we assume that the number of 
atoms in each bin is measured with some uncertainty, and convolute the probabilities 
entering the fit likelihood function with a Gaussian error distribution, p e rr{nk\n' k ) = 
e~( nk ~ n '^ 2 /( 2a c™/ m \ As a result, the variance is increased by a 2 VJ ./m. Very 

promising detection techniques based on the detection of fluorescence photons rely upon 
detection of on average a fluorescence photons per atom, and this number fluctuates at 
the shot noise level. Error propagation gives a 2 rr = -(rik) and the phase uncertainty (J7J) 
becomes 

a22M _ A 2 # (m) + — - - 

flt flt mav^Mrf' 

2^k=l (n k ) 

Using = 100 and the optimal state denoted by the vertical line in Fig. |2j for the 
bin size Ax = — a SSN phase uncertainty is preserved for a > 2.2 - a condition well 
satisfied by the light-sheet technique, which can give a ~ 10 photons per atom [25J. 
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Figure 4. The phase uncertainty (|13p. calculated for N = 100 in the optimal state 
(denoted by the vertical dashed line in Fig. |2}, as a function of a for four different 
values of 77. The horizontal dashed line denotes the SNL. The inset shows the phase 
uncertainty Eq. (fl2|) as a function of N (solid black line) with the optimal state, 
compared to the scaling from the Gaussian approximation, mA 2 0»Ji pt = 27V~ 4 / 3 
(dashed red line), the scaling of the QFI mA 2 ^[ = iV~ 4 / 3 (dotted blue line) and 

the scaling of Eq. (T3J} with a = ^ and 77 = 90%, giving mA 2 6^' opt = 1.48JV -1 - 16 
(dash-dotted green line). 



5. Conclusions 

We have derived an expression for the phase estimation uncertainty for a generic 
situation where the phase is inferred from the positions of probe particles, when only the 
one-body density is known. In a Mach-Zehnder type interferometer, the sensitivity of 
this protocol coincides with the well known error propagation formula and we recover the 
known Heisenberg limited phase uncertainty, A9 oc N~ l . Then we consider the simplest 
"double-slit" interferometer based on spatially interfering wave-packets suggested in 
Ref. [28], which still scales at best as A9 oc iV~ 2//3 , limited by the loss of fringe 
visibility. Nevertheless, the phase uncertainty for the interference pattern (|12l) has a 
major advantage over the MZI Q. Namely, it performs equally well for any value of 9, 
while ([9]) can reach very high values around 9 = | and 9 = |7T. 

The interferometric protocol employing the interference pattern could be 
implemented with a BEC trapped in a double-well potential. After imprinting the 
phase and switching the trap off, the two clouds would expand and interfere. The 
atoms could then be detected using for instance the light-sheet method [25], based on 
the fluorescence measurement of photons scattered by the atoms crossing a laser beam. 
Another possible scheme relies upon letting the atoms fall onto an optical lattice. If 
the interference pattern is dilute so that there is no more than one atom per site, their 
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positions could be detected by a fluorescence measurement [25] with ultra-high efficiency 
and resolution. Recent quantum interferometry experiments [91 [111 EH] indicate that 
very important effects limitating the phase uncertainty are the detection imperfections, 
which we believe to have realistically taken into account in our proposal. Another 
relevant constraint to the precision of a double- well interferometer comes from the noise 
present in the interferometric sequence i). A recent theoretical work |36j shows that 
the amount of squeezing at the optimal point = 0.44 for N = 100 particles) could 
be reached with a double-well BEC using a refocusing method even in presence of the 
latter source of noise. 
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Appendix A. Derivation of the phase uncertainty 

Appendix A. 1. Single-atom detection 

The proof that the estimator taken as the phase 9y^ which maximizes the likelihood 
function of Eq. ([2]) is unbiased and has the variance of Eq. ([3]) asymptotically in m 
can be performed along the lines of the standard proof that the Maximum-Likelihood 
method saturates the Cramer- Rao bound 1 39 J for m ^> 1 




where F is the Fisher information to be defined below. The original proof was given in 
Ref. [39]. Recent formulation can be found in Ref. [40], and a more accessible albeit less 
rigorous version is given in Ref. [4~T] . 

The proof consists of two steps. Firstly, we show that the estimator is consistent, 
i.e., for m — > oo, the probability that 9y^ ^ 9 goes to zero, which means that 
the estimator approaches the true value of the phase shift asymptotically. This 
demonstration follows Ref. [40]. Secondly, we adapt the simplified proof of Ref. [41] 
(see also [22]) that a consistent ML estimator is also efficient, which means that it 
saturates the Cramer- Rao bound, Eq. ([A. II) . for m ^> 1. In particular, the proof shows 
that for large m, 9^ is distributed with a Gaussian distribution with the variance of 
Eq. (2) and mean 9, hence it is also unbiased. 

Before we start, a remark is in order. It might seem surprising or even wrong to 
use the likelihood function of Eq. fl2]). The reason is that we use only the single particle 




A 9 




1 



(A.l) 
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probabilities to estimate the phase, even though in a single shot of the experiment, the 
N particles will be correlated in general. Hence the probability that in a single shot the 
second particle arrives at a position x<i generally depends on the position x% where the 
first particle was detected. In traditional ML estimation, one would therefore define the 
likelihood function as 

m 

£(^=n^(4V) (a.2) 

with the iV-particle conditional probability density pjv(x^|<^) instead and define 
the estimator as the maximum of this function. As mentioned above it can be 
shown [39], HOI HI] that this estimator is consistent, unbiased, and that it saturates 
the Cramer-Rao bound with the Fisher information 

F N = I dx ,1 , n J d ePN (x N \9)) 2 . (A.3) 



Pn(x n \6) 

This bound cannot be overcome by any other estimator using the results of the 
measurement governed by the probability density pn(xn\9) or by any of its reductions. 

However, there are no restrictions on how an estimator can be defined, there will 
only be differences in the performance of different estimators. We define the estimator 
based on the single-particle probability density only. As the proof below shows, this 
estimator is consistent, unbiased, and has the variance of Eq. (2). The variance is 
ultimately limited (but is generally larger than) the ultimate limit from Eq. (lA.lj) with 



the Fisher information from Eq. ( 1A.3I) . However, it has the advantage that it is accessible 
experimentally and also allows for sub shot-noise phase estimation, as shown in the main 
article. 

Consistency. We recall the definition of 9y^ , which is the value of the parameter <p 
which maximizes £{<p) from Eq. (1) from the main article. Equivalently, it maximizes 

1 m N 

/ (ra ^) = -£X>Pi(*?V), (A.4) 



m 

i=i k=i 



where N is the number of particles and m the number of independent repetitions of 
the experiment. The events xn are distributed with the conditional probability density 
Pn(xn\0), where 9 is the true value of the phase shift, and p\(xk\9) is obtained from 
Pn(xn\0) by integrating over all Xj^k- 

We assume identifiability, i.e., that pi(x\8) = pi(x\9') for all x is equivalent to 9 = 9'. 
Consistency is then proved by showing that f(<p) = lim /^(v?) has a maximum at 
<p = 9 as follows: 

1 m N 

m-f(o)= hm -j2J2( ln p^\v) -^pii^m 

m->oo m * — ' z — ' V / 



=1 fe=l 



N „ 

^ / dx N pi{x k \9)(\npi(x k \ip) -\npi(x k \9) 
fc=i J 

N [ dx Pl (x\9)\n^P^-<0, (A.5) 
J Pi{x\9) 
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where we have used ln(y) < y — 1. The equality is obtained iff y = 1. Hence the 
inequality Eq. flA.5j) is saturated iff pi(x\ip) = pi(x\6) for all x. It follows that tp = 9 by 
the identifiability assumption. Hence 9^ — > 9 for m — > oo. 

Efficiency. We expand the first derivative of £ from Eq. (1) of the main article 
around 9, 

d v log C{<p) ~ 9„ log £(</?) | fl + log £(</?) - 9). (A.6) 

We now set ip = 9y^. Since this phase maximizes £, the left-hand-side vanishes and we 
obtain 

,H m _ ^log£(<^)| e 



^log£(^)| e 



(A.7) 



The consistency of the estimator ensures that we can neglect terms of higher order in 
®ml ~ & provided that m is large enough. 

In order to investigate how 9^ —9 is distributed, we start by computing the average 
of the denominator, 

m N 



m 



i=l k=l 



N 



-►r»>l m E / ^NPN(&N\0)f%k>S(pi(?k\ ( P)) 

k=l J ' 8 

= —mNF x . (A.8) 

The coefficient N results from the indistinguishability of the particles and Fx is the 
Fisher information calculated with the single particle probability density, 

1 ,2 



Ft 



dx\- 



Coming back to Eq. I I A. 71) . we get 



(A.9) 



N 



□(m) 
ML 



(A.IO) 



1=1 1 1 fc=i 

Hence the difference Oy^ — 9 is the average of m random variables, which in the central 
limit are distributed with a Gaussian probability. With a calculation similar to the one 
of Eq. (1A.8j) one obtains that the average value vanishes, which means that the MLE is 
unbiased. The variance of the distribution in the central limit is 



A 9 



2flH 



ML 



mNF 1 



1 + (#-!)£. 

#1 



where 



C 



dx 1 dx 2 p 2 (xi,x 2 \9) 



d e pi{xi\9) d e pi{x 2 \9) 



Pi{x x \9) pi(x 2 \9) 
Therefore the Fisher information from Eq. (1A. 1|) is 

NF X 

1 + (N-1)% 



(All) 



(A.12) 



(A.13) 
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Hence the correlations enter via the two-particle correlation function p 2 even though in 
the definition of the estimator only the single body density p\ is used. 

Appendix A. 2. Multiple-atom detection 

The multiple-atom detection relies upon dividing the space into n^n bins, each of size 
Ax. In every bin, the number of atoms is measured and this result is averaged over m 
realizations, giving 

m 

n k = -Y^nf. (A. 14) 

According to the central limit theorem, for large m the probability of detecting is 
Gaussian, 

1 _(™fc-(™fc>) 2 

pfafcM = - jz A9 e . (A.15) 

Here, (n^) = lim fik and A 2 nfc are the associated fluctuations and both depend on the 

m— >co 

value of if. We construct the likelihood function as follows 

"bin 

C & t(tp) = Y[p(n k \<p). (A.16) 

fc=i 

Again, we follow the steps of the proof of Fisher's theorem and expand the derivative 
of the logarithm of the likelihood function around the true value and set ip = 9 ML , 



/ Mm) n\ 

For the following calculations we introduce the more compact notation 



(m) 
L > 

(A.17) 



d v u((p)\ e = d e u, (A. 18) 

used also in the main text. We start with calculating the denominator of Eq. ( 1A.17|) 
which in the large m limit reads, 

9 e log£ fit = -m g [-^ + [K> - nj + j 



K - (nk)] d 9 - 



(A.19) 



2 1 K x y A 2 n fc _ 

According to (IA. 18|) . both (n^) and A 2 nfc in the above equation are a function of 9. 
In analogy to Eq. (1A.8|) for m — > oo above denominator is replaced with its average 



value. Upon averaging, the second term vanishes and the third term is proportional to 
A 2 nk/m, which is negligible in the limit of large m. Therefore we obtain 

d^logU^-rn^^l. (A.20) 
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The average of the square of the nominator of Eq. (1A.17|) reads 

<(^log£ fit ) 2 > . m 2 £ 3^>|M (fo - - (n ,») 

^bin r\ I \ 

- ta, £^*^-w)(*-w)') 

^bin -i -i 

w E * " < " t>)2( ' 1 ' " ( " ,))2 )- (A - 21) 

k,i=i K 1 

The second term vanishes and the last, after the average is calculated in the central 
limit, becomes m-independent, and thus can be dropped when compared to the first 
term. 

Let us now separate the sum over k, I into k = I and k ^ I parts. The first one 
simply gives 

^ ( d e (n k )) 2 
A 2 n k ' 



k=l 

The non-diagonal part k 7^ I depends on the two-site correlations af^ = m(^n k — 

(n k ))(fii — (ni))\. To justify Eq. (5) of the main text, we notice that the atom- number 
fluctuations read 

A 2 n k = Pl (x k \9)Ax + (p 2 (x k ,x k \6) - Pl (x k \9) 2 )(Ax) 2 (A.23) 

and thus in a small bin-size limit are Poissonian, i.e. A 2 n k = P i(x k \9)Ax = (n k ). 
Therefore, 

"bin [de{n k )] 2 1 v-m bin „2 dg(n k ) dg{rn) 



E«bin [Qe{n k )\ , V^ n bin 
- = fc=l (n fc ) + 2^1=1 "fc,Z (n fc ) (n,> ^ ^ 



at 



ni 



Enun [de{n k )] 2 
k=l 



("fc> 

The last step is to calculate the cross-corelation term 

alj = (Ax) 2 [^(x^ixkie^ix^ix^)) - (&{x k \6)$(x k \e))(&(x l \6)${x l \e)) 

= (Ax) 2 N [(N - l)p 2 (x k , Xl \e) - N Pl (x k \9) Pl (x k \9)] , 
where in the last line we used 

(&(x k \e)*(x k \6)&(x l \6)i!(x l \6)) = 

(¥{x k \e)¥(x l \e)i>{x l \e)i>{x k \e)) = n{n -^{x^x^e), (A.25) 

which is true for k ^ I. 

Appendix B. Derivation of the phase uncertainty for the MZI 

When the two wave-packets are fully separated, so that ip a (x)ipb(x) = for all x, the 
one-body probability reads 

Pi(z|0) = ^(if\x\9)if(x\e)) = l - [\%{x)\ 2 {l -i/sin0) + |V6(^)| 2 (1 + ^sin 
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where v = jr(J x ) is the fringe visibility. This probability gives 



\Mx)\* + \Mx)\ A 



\^ a (x)\ 2 (l - usin9) + |^ 6 (x)| 2 (l + z/sini 



v 2 cos 2 6 
1 — v 2 sin 2 ^ 

Now we calculate the second order probability and get [cf. Eq. (jlj)] 

N 2 

N(N - l)p 2 {x 1: x 2 \9) = {^{x l \9)^{x 2 \9)^{x 2 \9)^{x 1 \9)) = — {\^ a {x x )\ 2 + \^ b {x x )\ 2 )p x {x 2 ) 



+ (|^ a (x 1 )| 2 " |V>6(Zl)| 2 )(|Va(^)| 2 - \MX2)\ 2 ) 



cos^((J 2 2 )-^)+sin^((J 2 )-^ 



+(|^(x 1 )| 2 - \M^)?Wa{^)\ 2 + \Mx2)\ 2 )j sm9(J x ). 

We insert this function together with p x (x\9) into the definition of C and obtain 

N(N - 1) • C = ^-(a + (3) 2 (J,)(/3 2 - a 2 ) sin9 

N 



+ ( (J 2 ) cos 2 6 + (J 2 ) sin 2 9 - j ) (a - /3) 2 



where 



a 



(J x ) cos# 



f -<J X ) sinfl 



and /3 



( j x ) cos 9 
f + (J x )sinfl' 



(B.l) 



If we now combine expressions for F\ and C as in Eq. ([3]), we obtain Eq. fl9]). 

Appendix C. The one- and two-body probabilites for the interference 
pattern 

The explicit expression for the one-body probability comes directly from the definitons 
of p\{x\9) and the angular momentum operators, 



Pl (a;|0) = —(&(x\0)if(x 



x 



b 2 



2 1 2(J X ) 

1 H jj— cos(kx + 



(C.l) 



Clearly, the coefficent v 



2{J X ) 
N 



is the visibility of the interference fringes. The two body 



probability reads [cf. Eq. (@|] 

(^{x 1 \9)^(x 2 \9)if(x 2 \9)^{x 1 \9)) 



p 2 (x x ,x 2 \9) 



N(N-l) 



1 - cos(k(xi - x 2 )) + [cos(/tXi + 0) + cos(kx 2 + 6)] 

4(J 2 ) 4(^ 2 ) 
+ ; ; ; ; ; - _ N cos(kx x + 9) cos(nx 2 + 9) + - — - sinf/txi + 6*) sin(/tX2 + 



N(N-l) 



N(N-l) 



(C.2) 
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Appendix D. Gaussian scaling 

To find the best possible scaling of the phase uncertainty ( 1121) with N, we model [V'in) 
with a Gaussian state as follows 

|<Ucxe-^^e~ ^ |j,iV-j>. (D.l) 

3 

The operator e _i ^ Ja: represents a beam-splitter which transforms the number- to the 
phase-squeezed state. For £^ = 1 the resulting state is spin- coherent, and by decreasing 
we increase the amount of phase-squeezing. 
This state is used to calculate the expectation values (J x ) and (J™) from Eq. f[T2"j) . 
For N ^> 1 the summation over j can be approximated with an integral. This way we 
obtain analytical expressions 

/V 2 3_ A /V 2 / 2 N 

{J x f = —e-^ and (J y 2 ) = — (l - e"^) . (D.2) 

which are then substituted into Eq. (TI2"]) . Taking ^ = , where < /3 < 1, we get 

mA 2 ^ = iV-^ +1 ) + N^r. (D.3) 

The phase uncertainty is optimal, when these two terms are equal, otherwise one of them 
would dominate at large N. This condition gives /3 opt = | and mA 2 ^ L '' opt = 2iV~ 4//3 . 
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